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ABSTRACT 

We report the discovery of a strong emission hne near 24.8 A (0.5 keV) in 
the newly discovered X-ray binary system MAXI J0556-332 with the reflection 
grating spectrometer onboard the XMM-Newton observatory. The X-ray light 
curve morphology during these observations is complex and shows occasional 
dipping behavior. Here we present time- and rate-selected spectra from the RGS 
and show that this strong emission line is unambiguously present in all the XMM 
observations. The measured line center is consistent with the Ly-a transition of 
Nvil in the rest frame. While the spectra contain imprints of absorption lines 
and edges, there appear to be no other significantly prominent narrow line due to 
the source itself, thus making the identification of the 24.8 A line uncertain. We 
discuss possible physical scenarios, including a gravitationally redshifted O vill 
Ly-a line originating at the surface of a neutron star or an unusual donor with an 
extremely high N/0 abundance (>57) relative to solar, that may have produced 
this comparatively strong emission line. 

Subject headings: accretion, accretion disks — binaries: general — X-rays: bina- 
ries — X-rays: individual (MAXI J0556-332) 



Introduction 



Monitor of All-sky X-ray Image (MA XI) detected a new X-ra y source on 2011 January 
11, subsequently named MAXI J0556-332 ( iMatsumura et al.ll201ll . hereafter we will refer to 
the source as J0556). The discovery was immediately followed up by pointed observations 
in other wavelengths. X-ray observations made with Swift indicated that the interstellar 
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extinction toward J0556 is small (Nh~10 cm , iKennea et al.ll201ll ) compared to most 



Galactic X-ray binaries (XRBs), which typically have Nh~10 cm~ . The low Nh, largely 
owing to the high Galactic latitude of the source (6=-25.183), implies that minimal line-of- 
sight absorption will facilitate efforts to study its low-energy X-ray spectrum. The X-ray 
spectral and timing properties of J0556 suggest that it is an XRB currently going through an 
episode of active accretion. The optical counterpart for this source was found to be a stellar 
object with USNO Bl.O ma gnitudes of m/j=19.91 and 771^=19.52, and an optical spectrum 
obtained by iHalperru (120111 ) showed emission lines from Ha, He I, and He H, indicating the 
presence of an accretion disk near the compact accretor. Spaceborne X-ray and groundbased 
I-band monitoring (using the CTIO/SMARTS 1.3m telescope) shows that the source is still 
active as of 2011 October 20. 

The nature of the compact accretor or the companion donor star in J0556 is currently 
unknown, although most evidence point towards a neutron star (NS) accretor. No "type 1" 
bursts have been seen from the source so far, but the X-r ay timing properties and evolution 
on color-color diagram are similar to those of NS XRBs (IHoman et al.l 120111 ). The optical 
to- X-ray and radio-to- X -ray flux ratios are also similar to Galactic NS XRBs (IRussell et al. 
201ll : ICoriat et al.ll201ll ). The distance to the source is also currently unknown. Analyzing 
the archival data of the source from the Catalina Real-Tir ne Transient Survey r evealed that 
no optical outburst had been detected since August 2005 (IMahabal et al.l 1201 ll ). Data from 
RXTE/All Sky Monitor show that this is its first X-ray outburst since January 1996. 



2. Data Reduction and Analysis 

XMM-Newton observed J0556 on MJD 55608 and 55653 (hereafter referred to as Obsl 
and Obs2 respectively), each observation being of about 40 k sec duration. We analyz ed the 



data obtained by the reflection grating spectrometer (RGS; Iden Herder et al.ll200ll ) using 
XMM-Newton Science Analysis Software (SAS; v. 11. 0.0) and following the standard extrac- 
tion procedures outlined in the XMM-Newton User Guid^. The background light curves 
were examined for the presence of proton flares. While there were no flares during Obsl, 
data from the initial 14 ksec of Obs2 had to be discarded due to high background. The SAS 
extraction regions for source and background spectra were inspected using the spatial- and 
energy-dispersion plots and left unchanged from the SAS default. The observed fluxes in 
each RGS CCD were compared with their 2% pile-up flux limits (§3.4.4.8.1 of the XMM- 
Newton User's Handbook), and we found that pile- up was significantly less than 2% even for 
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the steady- high state (when the source was brightest among the two observations). 

The average RGS count-rates during Obsl and Obs2 were 41.2 counts/s and 12.1 
counts/s respectively. The hght curve morphology during Obsl (Fig. [T]) is quite complex, 
showing dipping behavior for ~20 ksec during the middle of the observation as well as ex- 
tended epochs without any dips. This confirms previous sho rter X-ray observations of J0556 



wher e the dips were seen occasionally (e.g. as reported by IStrohmayed 1201 ll : iBelloni et al. 



20111 ). Given the marked changes in the source behavior, it is likely that the spectral prop- 
erties were changing during Obsl. Therefore we extracted spectra from 3 separate intervals 
based on the combined RGS count-rate [cps) and time [t seconds since the start of the ob- 
servation): (1) steady-high - when l,000<t<ll,000 and cps>45, which constitutes most of 
the time when the source flux was steady; (2) unsteady-high - when t>12,000 and cps>42, 
selecting the photons collected during the peaks of the dipping period, and (3) unsteady- 
low - when t>12,000 and cps>35, selecting the photons collected during the minima of the 
dipping period. These selections are marked in Fig. [1] by differently colored regions. Once 
the background flare was excluded, the Obs2 light curve was steady and did not show any 
dips or any other signs of large variability. Therefore we created one single time-averaged 
spectrum of the entire Obs2 after excluding the initial flare. 

We analyzed data from the full RGS bandpass (7-36A) for both RGSl and RGS2 units. 
Since the second order RGS spectra do not cover the wavelength region near the 24.8 A line, 
we restrict our analysis to the first order spectra only. The spectra were grouped for all fits 
with the criterion that there were at least 30 counts /bin. Fits to the X-ray spectra were 
performe d in waveleng th space using ISIS (v.1.6.1-36: iHouck fc Denicolal 120001 ) and XSPEC 



(vl2.7.0; lArnaudlll996l ). The solar abundance table given by lAnders fc Grevessd (119891 ) were 
used for all fits. 

Apart from the 24.8 A emission line which is readily visible in the RGS spectra (Fig. 
there are also many low-significance "bumps and wiggles". However there appear to be no 
strong narrow spectral features associated to the source that are of comparable strength 
to the 24.8 A line. Signatures of the intervening interstellar medium (ISM) are seen via 
the presence of neutral and ionized atomic oxygen edges and absorption lines in ~22. 1-23.8 
A region. Also evident in the brighter observati ons is the instrurnental neutral O K-edge 
absorption feature near 23.5 A as noted, e.g., by Iden Herder et al.l ( 200ll ). There is a hint 
of an excess near 29 A, which may be due to a blend of He-like Nvi resonance line (28.78 
A in rest frame) and intercombination lines (29.082 and 29.084 A). Note that the possible 
Nvi feature does not show the forbidden component at 29.535 A. But that is not surprising 
because the emission region in a compact XRB would be expected to be much denser than 
10^° cm^^. Also, if the excess is indeed due to Nvi, this would support a nitrogen-rich donor 
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scenario. While there are no other narrow hnes of comparable strength, there is a broad 
feature near 12.2 A (~1 keV; covered only by RGS2 and not by RGSl). The feature can be 
well modeled by a Gaussian centered near 12.2 A and a width of 1.4 A . Such a broad feature 



near ~1 keV is sometimes seen in XRBs (see, e.g.. ICackett et al.ll2010l ). Given the quality of 
the data, we cannot rule out the possibility that it could be instrumental, or relativistically 
broadened Ne X line, or smeared Fe L-edge seen in reflection. 

Various low-significance features are not only difficult to identify but also make the 
determination of the continuum difficult. For all the observations, the RGS continuum can 
be well modeled with a phenomenological broken power law model. We also attempted to 
model the continuum with a thermal (both blackbody as well as accretion disk models were 
tried) plus power law model. However the thermal component in these models was not 
constrained and the best-fit temperatures were much higher than that of the energy range 
considered. Thus in effect the thermal-|-power law models were mimicking a broken power 
law. The exact choice of continuum in the comparatively narrower RGS energy range is not 
crucial for measuring the 24.8 A line parameters, and we stress that the choice of the broken 
power law model is nothing but an empirical description of the continuum for the purpose of 
measuring the line. A full analysis of the 0.6-10 keV EPIC data is beyond the scope and focus 
of the current Letter, and will be presented in a later work (Maitra et al., in prep.). However 
we would briefly like to note that the phenomenological thermal-|-power law flts to the PN 
data shows no strong evidence for hardening of the power law component during the dips. 
The power law photon index during dips (unsteady-low) is r=2.65±0.10 and that outside 
dips (unsteady-high) is 2.51±0.06. However, both the temperature and normalization of the 
thermal component are smaller during the dips (when kT55=0.70±0.15 keV, Nbb=90±21) 
when compared to those outside the dips (kTBB= 0.81±0.02 keV, Nbb= 143±13). As a 
result, the ratio of unabsorbed non-thermal-to-thermal flux within the PN energy range is 
higher during the dips compared to that outside the dips. We modeled absorpt ion of X-rays 



by th e ISM using the Tuebingen-Boulder model (tbnew; an updated version of IWilms et al. 



2000l ). The best-flt model parameters for the various observations/selections are given in 
Tables [1] and [2l The quoted errors for the best-flt models correspond to a 90% confldence 
limit for the continuum model parameters, and 68% confidence limit for line parameters. 

While fitting, the normalization of RGS2 was allowed to vary w.r.t. that of RGSl. The 
fits show that the RGS2 normalization factor is in the range of 0.94-0.97. Fig. [3] shows the 
24.8 A line fiux and 26-28 A continuum fiux for the various observations/selections. Note 
that the ratio of line-to-continuum fiux during Obs2 (0.11±0.03) is similar to that during 
the steady-high state of Obsl (0.14±0.02). This strongly suggests that the line and the 
continuum fiuxes are causally linked and both originate from the accreting system (therefore 
ruling out the possibility that the line is from some other astronomical source in the angular 
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proximity of J0556; furthermore there is no other source of comparable X-ray brightness in 
the EPIC image). 



3. Abundance estimates using XSTAR photoionization code 

In order to obtain a self-consistent understanding of th e elemental abundances in the 



observed spectra we used the XSTAR photoionization code (iKallman &: Bautistdl200l[ ). A 
grid of 480 XSTAR models was created from the following parameters: (1) model column 
density of the emitting plasma was sampled between lo^^"^^ cm~^, (2) the ionization pa- 
rameter (log(^) where ^=L/nR^) was sampled between 0-4, (3) N abundance was sampled 
between 1-100 times relative to solar, and (4) O abundance sampled between 0.01-1 times 
relative to solar. A hydrogen nucleus density of 10^^ cm~^, and a covering fraction of 0.5 were 
assumed. A high turbulent velocity of 1,000 km s~^ was used to simulate the line broaden- 
ing. Phenomenological blackbody+power law fit to the 0.6-10 keV EPIC-PN data obtained 
during the Obsl/steady-high state (A;Tbb=0.86 keV, Abb=147, rpL=2.20, ApL=0.42) were 
used as user-defined in cident radiation field in XSTAR. For a source distance of 1 kpc (note 



that iWelsh et al.l 120101 reported a cloud of Na I and Ca II at ~100 pc in the direction of 
J0556; moreover, given the high galactic longitude of J0556, a distance of ^1 kpc would 
put the source far out of the galactic plane), the fiux measured by PN implies a luminosity 
of 2.6x10'^^ erg s~^. The remaining XSTAR parameters were set to their default values. 
An XSPEC table model created from this grid was fit to Obsl/steady-high state. The best 
fit suggests a model column density of 3.2x10^^ cm~^ and log(^)=2.3. The 90% confidence 
lower limit on log(^) is 2.1 (mainly based on an upper limit on the Nvi/Nvil ratio), but 
the upper limit on log(.^) could not be constrained from the data. Since there is only one 
statistically significant narrow line in the spectrum, it is not possible to constrain the abun- 
dance of any single element. Rather, the fits are only able to give limits on the relative N/0 
ratio. The fits suggest N/0 overabundance >57 (>27 at 90% confidence) with respect to 
solar. This translates to an absolute N/0 of >7.5 ( >3.6 at 90% confidence) since we use the 



solar abundance table of I Anders fc Grevessd (119891 ). 



We also created XSTAR models assuming the line is due to redshifted O vill Ly-a and 
solar nitrogen abundance. The best-fit to this redshifted O Vlll model with x^/z/=4702/3608 
is statistically slightly less favored than the N-rich model with x^/z/=4673/3607. We note 
that the best-fit log(^) was pegged at the maximum value allowed by our redshifted oxygen 
table model, implying a rather high log((^)>6 would be needed in that scenario. 
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Discussion and Conclusions 



J0556 is one of the few exceptional XRBs whose low-energy X-ray spectrum can be 
studied in detail due to low interstellar extinction. Fits to its RGS spectra imply a column 
density of (2.1-4.6) x 10^° cm~^, which is at least an order of magnitude below most XRBs. 
Both XMM-Newton observations of J0556 show a strong emission line near 24.8 A. If the 
line is from Nvil, this would require the donor to have an extremely high N/0 abundance 
(>57) relative to solar, based on the weakness of oxygen lines of similar charge states. 
Typically, for an XRB with a solar- type donor, the oxygen abundance is ~8x greater than 
that of nitrogen, and the X-ray spectra show lines from H- and He-like ions of oxygen. 
Even in the intermediate- mass X-ray binary Her X-1, wherein N/0 is 4x solar, O lines are 
observed to be as strong or much stronger than N lines from similar charge states (see, e.g., 
Jimenez-Garate et al.ll2002l ). While N overabundance is not uncommon in low- mass XRBs 
(e.g., in 4 of 5 UV spectra analyzed by iRaymondl ( 119931 ) the N/C ratios were several times 
solar, with a 9:1 ratio in Cyg X-2 being highest), the extremely high N/0 in J0556 makes it 
an unique XRB. 

Given the lack of source distance and any signature of orbital periodicity so far, it is 
only possible to speculate on the nature of the donor star in this system based on its strong 
N/0 overabundance and color (see below), and exotic stars such as hot, core- helium burning 
subdwarfs (sdB, sdO), or degenerate white dwarfs (WD) app ear to be strong cand i dates . Hot 
subdwarfs were first discovered at high Galactic latitudes by lHumason fc ZwickyI (Il947l ). and 
subsequent studies have shown that about half of them are in binary systems with WD or 
low-mass main sequence stars. It is generally thought that the hydrogen-rich outer envelope 
of the progenitor hot subdwarf in a binary system is lost via either Roche-lobe overflow or 
common-envel ope ejection mec h anism, thus e x posiii g CNO products dredged up from the 
core (see, e.g., iHan et al.ll2002l ). iNaslim et al.l (120101 ) have recently determined abundances 
for 6 sdB stars. The absolute N/0 ratios were 10, 13 and 20 for three of them (SB 21, 
LB 1766, and BPS CS 22940-0009), and oxygen was not detected in the other three. These 
stars were Ne-rich as well, so future observations of the Ne spectrum of J0556 would be 
especially interesting. One particular sdB star, PG 1219+534, is observed to have N/O>100 



(see, e.g., iHeberl l2009l for a review). Based on recent evolutio n ary c alculations studying 
formation of ultracompact binaries (UCBs) by iNelemans et al.l (|2010[ ). it is also possible 
that the donor is a helium WD. In the absence of detection of any C lines, N/C ratio is 
unconstr ained but likely high g iven the strength of the Nvil line. Based on the diagnostics 
given by iNelemans et al.l ( l2010l ). a high N/C would also point towards a helium WD. 



Assuming that the magnitudes reported in the USNO-Bl.O catalog for the optical coun- 
terpart were made during quiescence, the temperature inferred from the de-reddened B-R 
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color is ~46,000 K. This almost certainly rules out typical donors found in low-mass XRBs, 
which tend to be of late spectral type (K and M) and hence cooler. The high temperature 
suggested by the USNO photometry would also support a hot subdwarf donor scenario. A 
caveat to the above scenario is that the mass transfer rate during quiescence may also be 
non-negligible, so that there exists an accretion disk even during X-ray quiescence (making 
the optical colors appear blue). Also, He-sdOs are thought to evolve from merger of WDs, 
and it may be difficult to place such a merger product in a binary with a NS. The low column 
density towards the source makes it unlikely that the donor is a Wolf-Rayet star because 
Wolf-Rayet stars drive massive winds creating high internal columns, e.g., as seen in Cyg 
X-3. 

If the donor is indeed a WD/sdB/sdO, then the system is likely an UCB with or- 
bital period ranging from a few tens of minutes to ~ few hours. Observations of such 
UCBs are sparse, and their evolution is only recently being explored thoroughly (see, e.g., 
Tauris fc vari den Heuvell l2006l : iNelemans et al.l 120101 ). From a sample of 50 sdB binaries 
Geier et al.l (120081 ) found that possibly four could have a neutron star or black hole donor. 
However, des pite a sensitive searc h, no radio pulsations have yet been detected from these 
four systems (ICoenen et al.l 120111 ). If the donor in J0556 is a hot subdwarf, this would be 
the first direct evidence of a binary system with a hot subdwarf donor and a compact (black 
hole or neutron star) accretor. 

Another intriguing possibility, assuming solar-abundance plasma, is that the observed 
line is a gravitationally redshifted O vill Ly-a line (rest-frame A=18.967 A) originating from 
the surface of the NS. A photon emitted with a rest-frame energy of Eq from the surface of 
a slowly rotating, spherically symmetric NS of mass Mns and radius -Rns is observed by an 
observer at infinity to have an energy E where E/Eq = a/1 — (2G'/c^)(Mns/-Rns)- Based on 
the most precise mass estimat es from pulsar timing e xperiments, observed NS masses range 



between 1.25-2M0 (see, e.g., Kramer fc Wexl l2009l : iDemorest et al.l 120101 ). Assuming the 



observed 24.8 A line in J0556 is gravitationally redshifted O Vlll Ly-a, the above mass range 
would imply a radius of 8.9-14.2 km. In this redshifted oxygen line scenario, the width of the 
line could be used to put an upper limit on the relative size of the emission region as follows: 
assume that the width of the observed line profile is mainly due to superposition of lines 
originating at various heights above the NS. Since a photon emitted from the surface (i.e. 
from a radius of Rns) is more redshifted than a photon emitted from Rns+A/?, the observed 
width of the narrow line may be related to AR as Ai?/i?Ns ~ o"/(-E'o — E) ^ 0.023. This 
is obviously an upper limit to AR, since it ignores any bulk motion. However, a redshifted 
O VIII line scenario would also require other redshifted lines of oxygen to be present in the 
spectra, which are not seen. Also, the hint of Nvi resonance and intercombination lines 
favors the N/0 overabundance scenario. Planned high-resolution optical/UV spectroscopy 
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will be key in understanding the nature of this unique source. 

We thank the anonymous referee for helpful comments. It is a pleasure to thank Norbert 
Schartel and the XMM-Newton planning team for carrying out the ToO observations. We 
would also like to acknowledge the use of the daily monitoring data obtained by the RXTE 
which has provided excellent all-sky coverage of X-ray sources over the past 15 years. 
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Table 1. Best fit continuum parameters for tfie RGS data 



Obs/Selcction 


NhxIO^o 
(cm^^) 




(keV) 


r2 


AT 


Obsl / steady-high 


■i-»»_0.74 


J-S'-0.28 


n 475+0.019 
'.'■^'3-0.012 


^•^'-'-0.05 




Obsl / unsteady-high 


2.29tO-J5 


1 1O+0.03 
-'-■-'-S-0.03 


1 401+0.050 


^■°'J-0.15 


OQO + 0.003 


Obsl / unsteady-low 


-0.65 


r, c:o+0.42 
^■^^-0.98 


306+0.022 


, c:7+0.05 
-^■O'-0.05 




Obs2 


q qq+0.73 

o.OO_Q 4g 


^■'jy-0.56 


n qqfi+0.050 


1 qi+0.05 

-■-■•J -"--0.04 


U.U4i_(j Q26 



^Errors are 90% eonfidence limits. Col. (1) gives the observation and time/rate selection 
as discussed in ^ col. (2) is the best-fit Njj in units of 10^" atoms cm^; cols. (3,5) are 
power law photon indices for E < Ei, and E > Ef, respectively; col. (4) is the break point 
for the energy; col. (6) is the power law normalization at 1 keV, i.e., photons keV~^ cm"-^ 
s-i at 1 keV.; 



Table 2. Best fit narrow line parameters for the RGS data^ 



Obs/Selection 


Center 
(A) 


a 

(A) 


lO^xFlux 
(photons cm~^ s^^) 


Equiv. width 
(mA) 




Obsl/ steady-high 


24.801«;«1? 


n ifi«+o.020 


Z.yD_o 23 


190li^ 


4650.2/3609 


Obsl/ unsteady- high 


24.8001H11 


n r,i 1+0.026 
u.ziI q Q22 


d. l'»_o 21 


2731}^ 


5138.0/3902 


Obsl/ unsteady-low 


24.824ini^ 


+0.026 
U.l (0_Q 023 


4 1Q+0.31 
*-J^''-0.30 


3431^5 


2632.1/2575 


Obs2 


24.738in^^ 


„ 210+0.053 
u.ziy_o 03g 


U-DD_o 07 




4213.0/3580 



^Errors are 68% confidence limits. Col. (1) same as in Table [T] cols. (2-5) best fit Gaussian line center, 
width (ct), line flux, and equivalent width for the narrow line; col. (6) x^/d-o.f for the best fit to the 7-36 A 
RGS data. 




Fig. 1. — RGS light curve of the first XMM-Newton observation (Obsl) showing its complex 
morphology, especially epochs of dipping and non-dipping behavior. The variously colored 
regions represent the intervals from which spectra were analysed. 
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Fig. 2.— XMM-Newton/RGS spectrum of J0556, in 7-36 A range, obtained on MJD 55608 
during the steady-high period (i.e. the cyan region in Fig. [T]). Top panel: RGSl data in 
red and the RGS2 data in blue (binned for visual clarity only). Data from both RGS units 
were fitted jointly as described in text. The joint best-fit models for RGSl and RGS2 
are shown by red and blue histograms respectively. Rest-frame energies of few different 
lines/edges typically prominent in this range are labeled. The energies of the resonance (r), 
intercombination (i), and forbidden (f) lines are also shown. The broad feature in the 10-14 
A range was modeled with a Gaussian centered at 12.2 A and width (a) of 1.4 A. Middle 
panel: Residuals in units of standard deviation, when the normalization of the narrow line 
at 24.8 A (0.5 keV) was set to zero. Bottom panel: Residuals in units of standard deviation, 
to the best-fit continuum + line model. 



- 14 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



5.0 



4.5 



4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



1.0 



0.5 



Obsl /steady-high i — • — i 

Obsl /unsteady-high i — e — i 

Obsl /unsteady-low i — ■ — i 

Obs2 I — B — I 



I I ' ' I ' ' ' ' I I I I 



10 



12 



14 



16 



18 



20 



22 



26-28 A continuum flux (10'^ Photons/s/cm^) 



Fig. 3. — Photon flux from the 24.8 A emission hne vs. 26-28 A continuum flux during the 
various observations/selections. 
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Fig. 4. — Supplementary figure showing the multi- wavelength monitoring of J0556. The 
top two panels show H- and I-band light curve obtained with the CTIO/SMARTS 1.3m 
telescope. The next three panels show the 1.5-12 keV RXTE/ASM light curve (third panel 
from top), the 2-20 keV MAXI light curve (fourth panel from top), and the 15-50 keV 
Swift/BAT light curve (bottom panel). 



